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We investigate the behavior of the recently proposed quantum Google algorithm, or quantum 
PageRank, in large complex networks. Applying the quantum algorithm to a part of the real 
World Wide Web, we find that the algorithm is able to univocally reveal the underlying scale- 
free topology of the network and to clearly identify and order the most relevant nodes (hubs) of 
the graph according to their importance in the network structure. Moreover, our results show 
that the quantum PageRank algorithm generically leads to changes in the hierarchy of nodes. In 
addition, as compared to its classical counterpart, the quantum algorithm is capable to clearly 
highlight the structure of secondary hubs of the network, and to partially resolve the degeneracy 
in importance of the low lying part of the list of rankings, which represents a typical shortcoming 
of the classical PageRank algorithm. Complementary to this study, our analysis shows that the 
algorithm is able to clearly distinguish scale-free networks from other widespread and important 
classes of complex networks, such as Erdos-Renyi networks and hierarchical graphs. We show that 
the ranking capabilities of the quantum PageRank algorithm are related to an increased stability 
with respect to a variation of the damping parameter a that appears in the Google algorithm, and to 
a more clearly pronounced power- law behavior in the distribution of importance among the nodes, 
as compared to the classical algorithm. Finally, we study to which extent the increased sensitivity of 
the quantum algorithm persists under coordinated attacks of the most important nodes in scale-free 
and Erdos-Renyi random graphs. 
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I. INTRODUCTION 

It is of great interest to explore and classify the large 
amount of information that is stored in huge complex 
networks like the World Wide Web (WWW). A central 
problem of bringing order to classical information stored 
in networks such as the WWW amounts to rank nodes 
containing such information according to their relevance. 
A highly successful and nowadays widespread tool for 
this purpose has been the PageRank algorithm fH 1^, 
which lies at the core of Google's ranking engine. In the 
foreseeable future where large-scale quantum networks 
have become a reality, classifying the quantum informa- 
tion stored therein will become a priority. It is in this 
sense that the recently introduced quantum PageRank 
algorithm [3] is an important achievement as it consti- 
tutes a quantization of the classical PageRank protocol. 
This new quantum algorithm has shown, applied to small 
networks, a striking behavior with respect to its classi- 
cal counterpart, such as producing a different hierarchy of 
nodes together, paired with a better performance. In this 
paper we investigate the properties of the quantum algo- 
rithm for networks which model large real-world complex 
systems. We also test the algorithm on real- world data 
stemming from a part of the WWW. 

Complex networks are more and more pervasive and 
essential in our everyday's life, and thus the importance 
of network science. Consequently, considerable research 
is devoted to analyze and understand networks like the 
World Wide Web, the Internet, networks associated to 
transportation and communication systems and even bio- 
logical and social networks. Starting with the seminal pa- 



pers by Watts and Strogatz on small-world networks [4] 
and by Barabasi and Albert on scale- free networks [5 , re- 
searchers realized that most relevant networks belong to 
a class known as small- world scale- free networks, as they 
exhibit both strong local clustering (nodes have many 
mutual neighbors) and a small average path length while 
sharing another important characteristic: the number of 
links of nodes usually obeys a power-law distribution (the 
network is scale- free). Moreover, it has been found that 
many real networks, including the WWW, are also self- 
similar, see [6 . Such properties can often be related to 
a modular and hierarchical structure and organization 
which is essential for their communication and dynam- 
ical processes [7^91 . On the other hand, this hierarchi- 
cal structure could explain the existence of nodes with a 
relatively large number of links (or hubs), which play a 
critical role in the information flow of the system. Hubs 
are also associated with a low average distance in the 
network. Several general reviews and books on complex 
networks are now available, to which we refer the reader 
interested in more information on this topic p!Qm2] . 

Networks considered in this paper are modeled by three 
classes of graphs: The first type are Erdos-Renyi random 
graphs [13 . These graphs are constructed by connecting 
a given set of nodes with directed edges, each one added 
according to a certain fixed probability. The second class 
of graphs are scale-free graphs which were introduced by 
Barabasi and Albert to model the WWW [H [14]. A 
graph is dynamically formed by a continuous addition 
of new vertices which are connected preferentially to ver- 
tices which already have a large degree. We consider here 
a version for directed graphs published in [15] . The third 
family are hierarchical graphs which are also scale-free 
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but their clustering and degree distributions are nega- 
tively correlated as hubs have a smaller clustering coeffi- 
cient than nodes with a lower degree. It has been found 
that hierarchical graphs constitute also a good model for 
the WWW [71 [81. 

In the next sections we apply the quantum Google 
algorithm proposed in ref. [3 to representatives of the 
three classes of graphs. Our focus lies on directed scale- 
free networks, and also on hierarchical graphs as they 
are good models for the WWW, but we in addition also 
consider Erdos -Renyi random networks as a reference in 
order to contrast the results which we find for the former 
models. 

In particular, two main fundamental questions will be 
addressed: 

1/ Does the quantum PageRank algorithm preserve 
the structure of a scale- free network? In other words, 
we ask whether the ranking distribution obtained by the 
quantum algorithm also follows the same pattern of node- 
importance as the underlying scale-free network, and 
does not get mixed up with the distribution correspond- 
ing to a random Erdos-Renyi network. Similarly, for a 
Erdos-Renyi network we study whether the behavior of 
the quantum Page Rank algorithm is intrinsic to these 
networks and to which extent it differs from the one for 
scale-free networks. 

2/ Can one improve on the information gained on some 
of the properties of a scale-free network as far as the 
quantum Google algorithm is concerned? Specifically, 
scale-free networks are robust against random distur- 
bances which are unavoidable in a noisy environment. 
However, scale-free networks are vulnerable to coordi- 
nated attacks. 

These questions are outlined in more detail through- 
out the paper and clear answers are given. We hereby 
summarize briefly some of our main results: 

i/ We find that the quantum PageRank algorithm de- 
tects each of the representatives classes of complex net- 
works, in the sense that the classification results show 
that the ranking of the nodes follow a similar law than 
the degrees of the network. In particular, for a scale-free 
network this distribution remains scale- free under the ap- 
plication of the quantum PageRank algorithm. This im- 
plies that the quantum PageRank algorithm is expected 
to be robust with respect to random external noise. 

ii/ Remarkably enough, the detection of the hubs for a 
network in the class of scale-free networks with the QG 
algorithm is clearly enhanced with respect to the classical 
PageRank algorithm. In particular, the quantum PageR- 
ank algorithm is able to more clearly reveal the existence 
of secondary hubs in scale- free networks. 

iii/ We find that the quantum PageRank algorithm 
and the resulting rankings are more stable than the clas- 
sical PageRank protocol with respect to the variation of 
the damping parameter a that appears in the Google al- 
gorithm. The low dependence on this parameter, which 
inevitably had to be arbitrarily tuned in the classical al- 
gorithm, provides the quantum algorithm with a higher 



objectivity of the importance rankings. 

iv/ Our study shows that the quantum PageRank al- 
gorithm displays for scale-free graphs a power law scaling 
behavior of the importances of the nodes. Furthermore, 
this power law behavior is more favorable than the one 
the classical algorithm exhibits: Indeed, a smoother be- 
havior is connected to a more harmonious distribution 
of importance among the nodes, which enables the al- 
gorithm to better uncover the structure of hubs in the 
underlying scale- free graphs. 

v/ The enhanced sensitivity of the quantum PageRank 
algorithm to structural details of the networks is related 
to an increased sensitivity of the quantum algorithm un- 
der coordinated attacks of the most important nodes in 
scale- free networks. 



A. Conceptual Setting of the Quantum PageRank 
Algorithm 

Motivated by the fact that in a near-future scenario a 
certain class of quantum network will be operative [16]- 
21 , but not yet a scalable quantum computer, in [3] 
a class of quantum algorithms to rank the nodes in a 
quantum network was put forward. The algorithms in 
this class must be compatible with the classical one. In- 
deed, existing projects for large-scale quantum networks 
contemplate using the backbone of existing communica- 
tion networks, upgrading them to include the quantum 
hardware in order to store and manipulate quantum in- 
formation. In particular, the directed structure of the 
graph must be preserved, a feature which is crucial to 
measure a node's authority. An instance of this class of 
algorithms was explicitly constructed and analyzed for 
graphs of small size and it was found to serve as a valid 
quantum counterpart of Google's PageRank algorithm 
(for details see ref. [3 ). 

To perform such a quantum task^ it is important that 
the ranking algorithm must in principle incorporate some 
of the quantum properties of the network, like quantum 
fluctuations, and be objective. Indeed, the latter prop- 
erty earned much of Google's PageRank's success and 
was achieved embedding in the algorithm the random 
walk of a surfer exploring the WWW based on simple 
sensible rules. Within the same line of reasoning, we set 
up a quantum walk based algorithm that mimics the ex- 
ploration of nodes in a quantum network, in the setting 
when these are represented by states of a Hilbert space. 
The simple rules are encoded in the quantum dynamics, 
and in doing so the quantum nature of the networks and 
the information stored is properly taken into account. 
In the setting where a fully fledged large-scale quantum 
computer is not yet available, a key property of the quan- 
tum algorithm is that it be efficiently simulatable on a 
classical computer, that is, it must belong to the compu- 
tational complexity class P. Our algorithm is based on a 
single particle quantum walk and thus efficiently simu- 
latable. 
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Furthermore, the fact that the quantum algorithm con- 
tains a quantum walk at its heart paves the way to the 
analysis of its dynamics from a purely physical perspec- 
tive: In our work, we analyze extensively the localization 
properties of the quantum walk contained in the algo- 
rithm applied to several classes of networks (see sect. 
IV) showing its effects on ranking. To do so, we intro- 
duce the new concept of a quantum Inverse Participation 
Ratio (IPR). This is a generalization of the classical IPR, 
which is the main quantity that has been used extensively 
to probe localization properties in the study of classical 
random walks. 

Besides its application to future quantum networks, 
the Quantum PageRank algorithm can indeed also be re- 
garded as a valuable quantum tool, that can be efficiently 
run on a a classical computer, to perform a "classical" 
task, namely to rank nodes in a classical networks. In- 
deed, embedding nontrivially the network connectivity 
structure in the quantum dynamics, our protocol turns 
out to show several features that improve those present 
in the classical algorithm. 

Here we do not concentrate on possible quantum 
speedups and a detailed resource analysis for the quan- 
tum algorithm, which lies in the computational complex- 
ity P of efficiently simulatable algorithms. Instead, in 
this work we focus on the advantages that ranking nodes 
in classical networks using quantum algorithms displays, 
such as an increased resolution in the structure analysis 
of scale-free graphs or an increased stability with respect 
to the variation of the damping parameter. 

B. Operational Summary of the Quantum 
PageRank Algorithm 

In this section we briefly review the quantum PageR- 
ank algorithm from an operational point of view, whereas 
the reader can find the details of the construction in [3 . 
This quantum PageRank algorithm satisfies all the prop- 
erties of this class and represents a valid quantization of 
Google's PageRank algorithm. A step-by-step illustra- 
tion of the quantum algorithm is presented in Fig. [T] 
In Google PageRank's algorithm the ranking is per- 
formed setting up a random walk on the network. The 
walk uses as a transition matrix, known as the Google 
matrix G. The Google matrix is the weighted sum of 
two transition matrices. The first walk is driven by a 
modified connectivity matrix E where outgoing links to 
all other nodes have been added to every node that has 
no outgoing link. The second walk is a simple random 
hopping that connects every node to any other node. Ac- 
cordingly, the Google matrix associated to a given graph 
is defined as follows: 

G:=ai3+ii^l, (1) 

where 1 is a matrix with entries all set to 1 and N is 
the number of nodes. The parameter a is known as the 
damping parameter. 



The ranking of the nodes is performed measuring the 
probability to find the walker on each node when the sta- 
tionary distribution / has been reached i.e. when GI = I. 

In the quantum PageRank algorithm the idea is to set 
up a quantum walk on the nodes of the network and to 
perform the ranking of the nodes measuring the proba- 
bility of finding the quantum walker on such nodes. 
The quantum walk is a quantization of the Markov chain 
underlying the classical PageRank algorithm and is set up 
using Szegedy's method [22 which allows one to take into 
account explicitly the connectivity structure, and the di- 
rectedness of the network in particular. One can quantize 
a Markov chain on a A^-vertex graph that has as tran- 
sition matrix, the Google matrix G. This is performed 
by introducing a discrete-time quantum walk embedding 
the stochastic N x N matrix G on the same graph. 

The Hilbert space is the span of all vectors repre- 
senting the N X N (directed) edges of the graphs i.e. 
n = span{|z)i|j)2 ,withz, j e N x N} = C^(g)C^, where 
the order of the spaces in the tensor product is meaning- 
ful here because we are dealing with a directed graph. 
For each vertex let us define the quantum state vector, 

N 

IV'i):=E\/^l^')il^)2 (2) 

k=l 

that is a superposition of the quantum states representing 
the edges outgoing from the j^^ vertex. The weights are 
given by the (square root of the) Google matrix G. 

The idea of the quantum PageRank satisfying all the 
properties that define the class of possible ranking al- 
gorithms, is to set up a quantum walk starting from 
the initial state |?/^o) = XliLi whose dy- 

namics is governed by the quantum evolution opera- 
tor U := 5'(2n — 1) where S is the swap operator i.e. 

s = j:f,k=i urnim and n := Ef=i miu 

Although the Hilbert space is clearly N'^ dimensional it 
can be shown that the dynamics takes place in an in- 
variant subspace which is at most 2N dimensional. This 
allows one to numerically treat networks with a larger 
number of nodes. 

The rankings of the nodes in the quantum network is 
calculated using the instantaneous quantum PageRank: 

I,{Pi,t) = {^o\U^^'\i)2{i\U^'\4>o)- (3) 

Another quantity, called the average quantum PageRank 
can be defined: 

(J,(P,)):=^Xi W,t)- (4) 

t=0 

Whereas this quantity can be shown to converge for T 
sufficiently large, the Instantaneous Quantum PageRank 
does not converge in time [23] . 

In [3 nontrivial features of the quantum PageRank were 
uncovered, such as the instantaneous outperformance of 
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FIG. 1. (Color online) Schematic outline and summary of the quantum PageRank algorithm as proposed in [3]. a) The 
internet can be thought of as a set of pages (nodes of a graph) connected by directed hyperlinks (edges of the graph). The 
classical (quantum) PageRank algorithm can be regarded as a single walker performing a directed classical (quantum) random 
walk on the graph, b) The connectivity structure of the graph, as described by the connectivity matrix C, is of paramount 
importance to perform the ranking of the importance of pages both in the classical and the quantum case. In the classical 
(quantum) case the walker performs an incoherent (coherent) walk according to a combination of two hopping processes along 
the graph. The dynamics is governed by the "Google matrix" C, which describes the dynamics as a combination of (i) hopping 
according to a (patched) connectivity matrix E of the graph (parameter a, see main text and [3] for more details) and (ii) a 
fully random hopping processes (parameter 1 — a), where each node is connected to all other nodes of the graph, c) In the 
quantum PageRank algorithm the Hilbert space is spanned by the set of directed links between all pairs i and j of nodes of 
the graph, as tensor product states \i)i\j)2- The initial state iV^o) as well as the coherent discrete time evolution operator 
for the directed quantum walk (see main text) are determined by the Google matrix G. d) Quantum fluctuations can lead to a 
reversal of the order of importances of pages at certain instances of (discrete) time (so-called instantaneous outperformance) , 
as well as on the average over longer times (average outperformance). The latter effect is reflected by changes in the ordered 
list of pages (nodes) when the importance of pages according to the quantum PageRank algorithm is compared to its classical 
counterpart, the list of classical PageRank values. Red, blue and green curves in d) show the instantaneous quantum PageRank 
of nodes # 4, 5 and 7 of the seven- no de-graph shown in c), which was explored in Ref. and leads to the quantum PageRank 
and classical PageRank lists displayed in e). 



the algorithm or the violation of the hierarchy as pre- 
dicted by the classical algorithm. 

These properties motivate us to investigate the per- 
sistence of these novel effects on larger complex net- 
works. This paper is organized as follows: In section 
|TI|we present an analysis on graphs of the scale- free type 
for networks with hundreds of n odes, and in particular 
for a real- world network. In III) we have extended the 



analysis of the quantum PageRank algorithm on Erdos- 
Renyi and hierarchical graphs. From this we conclude 
that the behavior of the quantum PageRank is charac- 



teristic of each type of complex networks, and moreover, 
it is different from the classical PR algorithm: the clas- 
sical ranking changes by quantum fluctuations. The new 
quantum dynamics incorporated in the task of ranking 
also raises other important questions on the properties 
of the quantum walk embedded in the algorithm such as 
the localization phenomena of the walker on the network 
is analyzed in section |IV[ while the stability of the rank- 
ing with respect to noise (or damping parameter) is stud- 
ied in section [Yl In section |Vl| we address the question of 
whether the power law behavior displayed by the classical 
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PageRanks is preserved by the quantized algorithm and 
compute the scahng exponent. In section |VII| se inves- 



tigate a very practical situation and study how sensitive 
the quantum PageRank algorithm is under coordinated 



attacks in scale-free graphs. Section |VIII| is devoted to 
conclusions and outlook, and an appendix provides some 
technical details. 



II. QUANTUM PAGERANK ON SCALE-FREE 
NETWORKS 



We will analyze the predictions of the quantum PageR- 
ank algorithm on complex networks. We will focus on 
random scale-free networks because of their widespread 
appearance and relevance in real- world applications. In 
the next section we will also deal with the important 
cases of random (Erdos - Renyi) and hierarchical net- 
works for completeness in order to check whether the 
quantum PageRank algorithm preserves the characteris- 
tics of different classes of complex networks. Moreover, 
the study of random and hierarchical networks results to 
be useful to confront its features with scale-free networks 
and to draw interesting conclusions. 

Random scale-free graphs [24^ J5 are ubiquitous in 
nature. They appear as good models of the World 
Wide Web |26j, airline networks [27 or metabolic net- 
works [28l[29], just to name a few. These are networks 
that display a small fraction of hubs, i.e. nodes with 
a high connectivity, a property that follows from the 
degree distribution P{k) that shows a scale-free behav- 
ior, P{k) ~ . Scale- free networks exhibit intriguing 
properties which have been studied extensively, such as 
robustness against uncoordinated attacks [3QH32] . good 
navigability ^33H35) and controllability j36ti38] . 

One of the first models proposed to describe scale- 
free networks is the preferential attachment model O [14] . 
In this model links are preferentially formed to already 
highly connected nodes. A random directed scale-free 
model for the WWW was also introduced in [39 and a 
generalization appeared in JTb\. To produce the char- 
acteristic power-law degree distribution of degrees, the 
models consider two main mechanisms: growth and pref- 
erential attachment. A graph is dynamically formed by a 
continuous addition of new vertices and each new vertex 
is joined to several existing vertices selected proportion- 
ally to their in and out degrees. The generalized model al- 
lows also the introduction of directed edges between two 
already existing nodes. In this work we use graphs cre- 
ated with this model as implemented in NetworkX [40 . 

We next discuss results originating from the applica- 
tion of the quantum PageRank to complex networks of 
the scale-free type with hundreds of nodes ((9(100)). To 
this end, we have performed a numerical analysis using 
the quantum PageRank algorithm. We find that the algo- 
rithm clearly identifies that the networks are of scale-free 
type. The algorithm is able to point out the most impor- 
tant hubs. This is a task already well performed by the 
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FIG. 2. (Color online) Comparison of the quantum and 
classical PageRank on a real network originating from the 
hyperlink structure of www.epa.gov 41 . One can clearly see 
how the hubs in the classical algorithm tend to concentrate 
nearly all the importance. The insets show that the quantum 
algorithm is capable to lift the degeneracy of nodes in the the 
low part of the list (see text in sect. |ll| 



classical PageRank. However, the quantum PageRank 
algorithm has improved ranking capabilities, not concen- 
trating all the importance on these few nodes. Indeed, it 
is capable of unveiling the structure of the graph high- 
lighting also the secondary hubs (see figures |3],[4)|5]and 
captions therein). 

Furthermore, we find that the hierarchy as predicted 
by the classical PageRank is not preserved. This is a 
property already found in ^ for smaller networks. From 
our study we are able to clearly conclude that far from 
being an artefact of choosing small networks, this results 
to be a generic feature of the quantum PageRank algo- 
rithm. We also found that the quantum PageRank is 
able to lift the degeneracy of the nodes that have a lower 
importance. This feature can be seen clearly in figure 
|2] where we analyze a subgraph of the WWW obtained 
by exploring pages linking to www.epa.gov and available 
from Pajek [41 . 



III. QUANTUM PAGERANK ON 
ERDOS-RENYI NETWORKS AND 
HIERARCHICAL NETWORKS 

A. Erdos- Renyi Networks 

In this section we will briefly introduce the Erdos- 
Renyi class of random networks and analyze the perfor- 
mance of the quantum PageRank algorithm applied to 
them. 

This class of graphs was introduced by Paul Erdos and 
Alfred Renyi more than fifty years ago [13l |42j 03] , and 
is of particular importance in the context of Graph The- 
ory. There are different equivalent methods to describe 
this family. To allow an easy computer implementation 
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FIG. 3. (Color online) Scale-free graph with 32 nodes and 
the comparison of the importance of the nodes when evalu- 
ated with the quantum and classical Page Rank. The classical 
PageRank shows a very sharp concentration of importance on 
the three main hubs, nodes 0, 1 and 2. One can see from the 
comparison of the predictions of the two algorithms the rel- 
ative emergence of secondary hubs (nodes 6, 7 and 14) when 
the importance is calculated with the quantum PageRank (see 
text in sect. [TH 
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FIG. 4. (Color online) Scale- free graph with 64 nodes and 
the comparison of the importance of the nodes when evalu- 
ated with the quantum and classical PageRank. The classical 
PageRank shows a very sharp concentration of importance on 
the three main hubs, nodes 0, 1 and 2. The quantum PageR- 
ank algorithm is able to better distinguish and to highlight 
the secondary hubs (in this graph nodes: 12 and 25) whose 
importance rises with respect to the primary hubs, (see text 
in sect. [lT| 
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FIG. 5. (Color online) Scale-free graph with 128 nodes (see 
text in sect. [ll| and a comparison of the importance of the 
nodes when evaluated with the quantum and classical PageR- 
ank. The classical PageRank shows a very sharp concentra- 
tion of importance on the three main hubs, nodes 0, 1 and 2. 
One can see from the comparison of the predictions of the two 
algorithms the relative emergence of secondary hubs (nodes 
8, 9, 10, 20 and 21) when the importance is calculated with 
the quantum PageRank (see text in sect. |ll|. 



we use the following: a graph of order N can be con- 
structed by connecting N vertices randomly by adding 
edges with a given probability, which is independent from 
other edges. We use directed versions of the graphs cre- 
ated with NetworkX [40 . 

The graphs falling into this class follow a Poissonian 
degree distribution, i.e., P{k) ~ (/c)^exp {—{k))/k\ where 
(k) is the average degree. Thus, most nodes have a degree 
not far from the average and therefore the graphs do 
not display relevant hubs when applying the PageRank 
algorithm. 

We have performed a numerical study on graphs of the 
Erdos-Reny type and found that the quantum PageRank 
algorithm displays a sharp change in hierarchy (see fig- 
ure |6]b). We can draw two fundamental conclusions: i/ 
The behavior of the quantum PageRank algorithm is in- 
trinsic and characteristic of the Erdos-Renyi class of net- 
works, in particular, the no-hub-behavior is reflected at 
the quantum level, ii/ However, the quantum PageRank 
algorithm changes the ranking of the classical algorithm, 
in this sense outperforming the classical PageRank hier- 
archy. 
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FIG. 6. (Color online) Quantum and classical PageRank 
in random (Erdos-Renyi) networks. Subfigure a shows a pro- 
totypical example of a random network of 64 nodes clearly 
indicating the absence of hubs in this class, b) A comparison 
of the importance as obtained from the quantum and classi- 
cal PageRank applied to a 64 nodes random graph. The im- 
portances calculated using the quantum PageRank algorithm 
display a change in hierarchy (see text in sect. III). 



B. Hierarchical Networks 

Some relevant real-life networks which describe tech- 
nological and biological systems, such as the WWW, 
some electronic circuits and protein or metabolic net- 
works are usually scale- free but have also a modular 
structure [6l[44]. That is, they are composed of mod- 
ules that group different sets of nodes. These modules 
can be distinguished by the fact that nodes belonging 
to the same module are usually strongly connected. On 
the other side, modules are relatively weakly connected 
among them. Thus, even when the networks are scale- 
free, their hubs use to have a low clustering as they joint 
different modules. Several authors claim that a signa- 
ture for a hierarchical network is that, other than the 
small- world scale-free characteristics, the scaling of the 
clustering of the vertices of the graph with their degree 
follows d (X l/ki [7, 8 . 

Hierarchical network models usually are constructed 
from recursive rules. For example, we can start from a 
complete graph Kn and connect to a selected root node 



n — 1 replicas of K^. Next, n — 1 replicas of the new 
whole structure are added to this root. At this step the 
graph will have vertices. The process continues until 
we reach the desired graph order. There are many vari- 
ations for these hierarchical networks, depending on the 
initial graph, the introduction of extra edges among the 
different copies of the complete subgraphs, etc. However, 
once the starting graph is given, these networks do not 
have adjustable parameters and their main characteris- 
tics are fixed. 

In [45], Barabasi et al. introduced a simple hierarchi- 
cal family of networks and showed it had a small- world 
scale-free nature. The model was generalized in [7 and 
further studied in [46] . For our analysis we have designed 
a directed version based on these graphs, see figure 
In this case the starting point is a directed 3-cycle. 

Another interesting family of hierarchical directed 
graphs has been obtained by giving directions to the 
edges of the construction published in [471 HH] , see fig- 
ure [7^. The graphs are in this case small-world, self- 
similar, unclustered and outerplanar (a planar graph is 
called outerplanar if it has an embedding where all ver- 
tices lie on the boundary of the exterior face). However, 
they are not scale- free, but follow an exponential distri- 
bution. It has been shown that many algorithms which 
are NP-complete for general graphs perform polynomial 
in outerplanar graphs [49 . 

We have performed a numerical study on hierarchical 
networks using the quantum PageRank algorithm. We 
analyzed two families of graphs (see figure [t] for the con- 
struction) and we find that the hierarchy (similarly to 
what was found in ref. [3 for the binary tree) is pre- 
served by the average PageRanks. Interestingly, though, 
the quantum PageRank is able to highlight the connec- 
tivity structure of the nodes that belong to the same level 
in the hierarchical construction (see figures [8]and[9|. We 
observe that the difference in importance between nodes 
belonging to the same hierarchical level but with differ- 
ent local connectivity is amplified when calculated using 
the quantum PageRank algorithm. 



IV. LOCALIZATION-DELOCALIZATION 
TRANSITION 

We have seen in the previous sections that the quan- 
tum PageRank algorithm is able to distinguish between 
networks of Erdos-Renyi and scale-free type. In particu- 
lar, the quantum algorithm, as opposed to the classical 
one, is also able to highlight the secondary hubs in the 
case of scale- free networks. Furthermore, regarding the 
classical and quantum algorithms as walks on a network, 
a necessary condition that few nodes, the hubs, account 
for the majority of the importance is that the walker 
be localized, i.e. the number of nodes with a significant 
probability to find the walker is negligible with respect 
to the number of nodes in the network. 
In this section we will, therefore, study the localization 
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FIG. 7. (Color online) Construction of hierarchical net- 
works, a) The family of outerplanar directed hierarchical 
graphs. Note that the generation labeled by n has 2^^^ nodes. 
We consider graphs of the generations with n = 4, 5, 6. b) The 
family of directed hierarchical graphs. In this case the gener- 
ation labeled by n has 3^ nodes. We consider graphs of the 
generations with n = 2,3,4. 



properties of the walker on different networks, of Erdos- 
Renyi and scale-free type. We will briefly review the case 
of the walker in the classical PageRank algorithm and 
explore which phases does the quantum walker choose in 
the case of the quantum PageRank as a function of the 
structural properties of the underlying network. 
In order to measure the localization phenomenon we will 
employ the Inverse Participation Ratio (IPR). This con- 
cept was introduced in the context of condensed matter 
and, more specifically, to study the Anderson localiza- 
tion phenomenon in disordered systems (see for exam- 
ple [50]). It is particularly useful to study localization- 
delocalization transitions. 

The IPR ill the case of a classical walker, is defined 



N 



(5) 



where r > is an integer parameter which can be freely 
chosen and is fixed. 

Let us consider a classical walker on a network. We can 
have two extreme behaviors. The first one is that the 
walker is completely delocalized i.e. the probability of 
finding it on a site is uniform. Therefore, if we introduce 
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81 nodes 



FIG. 8. (Color online) Comparison of the quantum and clas- 
sical PageRank for the family of hierarchical graphs described 
in the text (see section III B and figure [tJ) for the construc- 
tion). We consider graphs with n = 2,3,4 (see subfigures a, 
b and c respectively). We find that the quantum PageRank 
preserves the hierarchy of the nodes but in addition it is able 
to highlight the connectivity structure of the nodes belonging 
to the same level. Indeed, in subfigure a, for example, the 
difference in importance between nodes 2 and 3 and 5 and 6 
is amplified when calculated using the quantum algorithm. 



a random variable X whose realizations are the sites of 
the lattice, we can write Pr(X = i) = 1/N ^ Mi. The 
other limiting case is that the walker is localized only on 
one site, i.e. the probability of finding the walker is a 
Kronecker delta, that is Pr(X = i) = Sij_, if the walker is 
localized on say, site j. The IPR (eq. ([5|) yields for the 
two limiting behaviors: 



1 

^l-2r 



if the walker is localized 
if the walker is delocalized. 



(6) 



Thus the IPR, displaying respectively, a power law or a 
constant behavior as a function of the number of nodes, 
is a legitimate witness of the localization of the walker. 
We can rewrite £^ci as: £^ci — N~^'^'^ . In order to study 
the localization-delocalization transition, it is useful to 
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FIG. 9. (Color online) Comparison of the quantum and 
classical PageRank for the family of hierarchical graph de- 
scribed in the text (see section III B and figure [T^l for the 
construction). We considered graphs with n = 4,5,6 (see 
subfigures a,b and c respectively). Also in this case we find 
that the quantum PageRank preserves the hierarchy of the 
nodes highlights moreover the connectivity structure of the 
nodes belonging to the same level. 



PageRank vector /, is given by: 
GI = I 



(8) 



and represents the stationary probability distribution of 
the walker on the network. 

In ref. [51 it was found that delocalization is absent for 
a very large range of values of a ranging from 0.4 to 1. 
This is quite natural because for the range of values of 
the damping factor a stated above, the second term, cor- 
responding to random hopping in the Google matrix (see 
eq. ([T])) is suppressed. This yields the localization effect. 
On the other side, for a close to 0, it was found that 
the walker is delocalized over the network. This is un- 
derstood easily, indeed for a close to only the second 
term in the Google matrix survives and the walker is only 
subject to random uniform hopping between any pair of 
nodes. This leads to a trivial phenomenon of delocaliza- 
tion of the walker over the whole network. 

Let us now focus on the localization phenomenon in 
the case of the quantum walk on a network according 
to the quantum PageRank protocol. In order to carry 
out the analysis we need to generalize the definition of 
the IPR given above reinterpreting the notion of prob- 
ability of finding a walker on a node when dealing with 
the quantum PageRank. In doing so, we will choose as a 
guiding principle the interpretation of the average quan- 
tum PageRank of a node (see eq. Q) as the probabil- 
ity of finding the quantum walker on a particular node. 
Therefore, we will employ the definition given in eq. Q 
that we will rewrite as: 



(2) 



(9) 



where M> is the (strong) measurement operator on the 
second space (indexing the nodes where the edges point 
to, see subsect. 



IB), i.e. Mf^ = |z)2(i|. 



We are now in a position to define the IPR ^ in the case 
of a quantum walk: 



introduce the normalized anomalous dimension A2r' 



N 



(10) 



T2r ■■= (2r - 1) + A: 



2r5 



(7) 



which interpolates between the two phases if the system 
undergoes a transition from a localized regime (where 

= 1 — 2r) to a delocalized one (where = 0). 
In [51 the localization-delocalization transition for a clas- 
sical walker performing a random walk (with transition 
matrix given in ([T]) ) was characterized studying its de- 
pendence on the damping parameter a where < a < 1. 
This study is important to understand at a deeper level 
the classical PageRank algorithm. Indeed, as we have 
anticipated, in the case of a scale-free graph observing 
localization over a broad range of values for a is indeed a 
necessary condition for the algorithm, quantum or clas- 
sical, to perform well the task of ranking the nodes. The 



Also in the case of a quantum walk we have ^ = N~^^'^ 
and it is evident that one can extract the localization 
phase of a walker from the scaling exponent of the IPR 
as a function of the number of nodes N. Indeed, from 
equation: 



log^^(l-2r-A2.)log7V 



(11) 



it is clear that the witness of the localization lies in the 
slope of the graph of the aforementioned log-log plot. 

We consider two kinds of networks, of the scale- free 
and of the Erdos-Renyi type. In order to study the local- 
ization phenomenon we generated networks with different 
numbers of nodes belonging to the two aforementioned 
classes. We then calculate the IPR (in the following we 
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will fix the parameter r = 1) in order to understand 
whether for ol = 0.85 the quantum walker was localized 
or delocalized on the network. 

We find that the IPR in the case of the class of scale-free 



networks does not vary appreciably (see fig. 10 ) signaling 
localization of the walker on the graph. Notice however 
that also in ref. [ST for these values of a a similar behav- 
ior was found. 

We have analyzed also the graphs in the Erdos-Renyi 
class performing the same steps as above. Our study 
shows that, albeit this being random graphs and that 
for this value of the damping parameter a the walk is 
strongly influenced by the topology of the network, the 
quantum walker is delocalized in this case. Indeed it can 



be seen from figure 10 that the behavior of the logarithm 
of the IPR is linear in the logarithm of the number of 
nodes. This is a clear witness of the delocalization phe- 
nomenon. 

This behavior is remarkable for two reasons. The first 
reason is that in the classical random walk case local- 
ization of the walker was found [51 . The second reason 
is that albeit randomness being expected to give rise to 
localization, in the case of graphs of Erdos-Renyi type, 
our study shows that the opposite is true. 
We conclude that the scale-free graph seems to favor a 
localization phase in both walks, random and quantum, 
respectively for the classical and quantumPageRank al- 
gorithm. This result in the latter case was obtained with 
a value of a = 0.85. This is consistent with a good rank- 
ing of nodes in a network. Indeed, in order to unveil the 
main hubs the random or quantum walk must be able 
to localize the walkers on few important nodes. Interest- 
ingly, instead, the Erdos Renyi graphs seem to prefer a 
delocalized phase albeit the networks being grown ran- 
domly. Both for the classical and the quantum PageRank 
one finds delocalization. This can be correlated with the 
absence of a small number of main hubs in this class of 
networks. 

As for the classical PageRank the localization- 
delocalization transition was characterized as a function 
of the damping parameter, it is important to understand 
how the quantum PageRank depends on the value of a. 
We will study this dependence and characterize it in the 
next section. 



V. STABILITY OF THE QUANTUM GOOGLE 
ALGORITHM WITH RESPECT TO THE NOISE 
PARAMETER 

In the previous section we analyzed the localization 
properties of the quantum PageRank algorithm. We 
studied this phenomenon having fixed the damping pa- 
rameter using the value a = 0.85. The natural question is 
to study how the quantum PageRank varies with respect 
to the variation of a. 

The stability of the quantum PageRank is an important 
issue to consider also because the damping parameter is 



arbitrarily tuned to a specific value. Indeed there is no a 
priori argument to fix the value of a and the value 0.85 
was originally chosen in the classical PageRank protocol 
to mimic the behavior of a surfer (or a random walker) 
that randomly hops to any other page once every seven 
times. Only a posteriori it turned out that this is indeed 
a sensible choice given that the network is small- world, 
being in fact a crucial ingredient for the PageRank algo- 
rithm to yield reasonable ranking results. In view of the 
ad-hoc choice of the precise of the parameter a, it is a 
very desirable property that the output of the algorithm 
be stable, i.e. the ranking vary slowly with respect to the 
variation of the damping parameter. A question that was 
addressed in the computer science community has been to 
quantify how susceptible to changes in this parameter the 
classical PageRank algorithm is. It was found that the 
effect of this parameter on ranking is large and that two 
rankings obtained by running the algorithm using differ- 
ent values of this parameter can be very different [51 . 
To tackle this problem in the quantum case we will make 
use of two quantities. The first one is related to the 
quantum fidelity (see e.g. [52l|53]) that provides a way to 
measure the distance between two quantum states. The 
second quantity that we will use is the classical fidelity. 
It is employed for the same task when dealing with prob- 
ability distributions. 

In ref. [51 the stability of the classical PageRank was 
studied as a function of the damping parameter. Since 
the PageRank vectors are classical probability distribu- 
tions one can measure the distance between two PageR- 
ank vectors, calculated using different values of the 
damping parameter, with the classical fidelity. The latter 
can be written as: 



f{a, a') = Y,Jl{Pj,a)I{Pj,a') 



(12) 



It was found that there is a plateau around (a, a') ^ 
(0.5, 0.5) and that the fidelity /(a, a' = 0.85) is rather 
flat around a ~ 0.85 thus implying that the classical 
PageRank is rather robust against perturbations (see 
ref. [51 ). 

The quantum fidelity is a quantity that measures the dis- 
tance of two quantum states: 



(13) 



which in the case of commuting density matrices reduces 
to the classical fidelity. 

Another valid measure of the distance of two quantum 
states is the trace distance: 



D{cr,p) = ^tr|p-cr| 



(14) 



where |r| denotes the square root of the (positive) oper- 
ator rV. 

The fidelity and the trace distance turn out to be equiva- 
lent measures of distance. Indeed if the fidelity of quan- 
tum states is near to one then their trace distance is close 
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FIG. 10. (Color online) The IPR (for r = 1) for networks of different classes, using the classical and the quantum walk, plotted 
versus the number of nodes in a log-log scale, a) The IPR using the quantum walk in the case of a scale-free graph and b) an 
Erdos-Renyi graph, c) The IPR using the classical walk in the case of an Erdos-Renyi graph. The value of a is set to 0.85 in all 
cases. In order to obtain the data we created using networkx 4 networks in the same class and with the same parameters having 
32, 64,128 and 256 nodes. We then calculated the IPR. In order to infer the phase of the walker we plotted log^ vs. logA^. A 
constant behavior signals that the walker is localized whether a monotonically decreasing behavior signal delocalization over 
the network. We can infer that the quantum walker is in a localized phase in the case of the scale-free network (see (a)). 
This behavior is in contrast with what is displayed in the case of the Erdos-Renyi graph where delocalization is found in both 
the quantum and classical case (see subfigures (b) and (c)). This result holds true also when one considers graphs with a 
significantly lower link to node ratio. In this case we have superimposed two lines that result from a linear fit. These have 
equation log^ = aqlogN + hq with aq — —0.8565 and hq — —0.2482 for the quantum PageRank (see subfigures (d)). In the 
case of the classical PageRank the equation is log^ = ad \ogN + bd with ad = —1.0932 and bd = 0.7125 (see subfigures (e)). 
These results are consistent with delocalization (confront text in sect. IV). 



to zero and viceversa [54] Therefore, either can be used 
for our purpose of measuring the stability of the quan- 
tum PageRanks and the choice is a mere matter of con- 
venience. To tackle the problem of stability in the case 
of the quantum PageRank we will use the trace distance 
and the classical fidelity. 

Let us rewrite explicitly the definition of the averaged 
quantum PageRank {Iq{Pi^a)) adding the dependence 
on a (that enters in the initial state and in the evolution 
operator of the walk) as: 



(15) 



where p]^{2t) = U^^ p]^ {0)Ul'^^ and for bookkeeping pur- 
poses it has been made explicit to which spaces the den- 
sity matrix refers to. 

Let us discuss how to apply the concept of trace dis- 
tance in our case. We might measure the instantaneous 
distance D{pa{2t), pa'{2t)), thus measuring the distance 
of the quantum states. However, it is more significant to 
take the time-average of the distance between the par- 
tial traces of the states, because we are interested in 



the quantum PageRank as an observable, rather than 
the state itself. Not being interested in the distance of 
states^ we can make use of a less refined measure of dis- 
tance (that is, one that can appreciate less the difference 
of two states). This can be written in the following form 
(see the appendix for the derivation): 



= maji\{Ig{Pi,t,a)) - {Ig{Pi,t,a'))\ 



(16) 



To summarize, we will study the stability making use 
of the classical fidelity and the quantity in (16), that is 



a simpler measure descending from the trace distance 
(which turns out to be equivalent to the quantum fi- 
delity). We perform an analysis of how the quantum 
PageRank varies with respect to the rankings when the 
value of a goes from 0.01 to 0.98. We analyze a scale 
free with 128 nodes that was generated with networkx, a 
python module. The results show clearly that the quan- 
tum PageRanks vary very little when the quantum walk 
underlying the quantum PageRank has a different damp- 
ing parameter a. Indeed it can be seen from figure pT]t hat 
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FIG. 11. (Color online) Measures of distance of quantum Page Ranks obtained with different values of the damping parameter. 
The network analyzed is a scale free network with 128 nodes (generated with NetworkX). The damping parameter varies ranging 
from 0.01 to 0.98. a) The fidelity obtained by applying the classical fidelity (see eq. (12)). One can see that there is a plateau 
for values of a around 0.8. b) The measure of distance obtained from the trace distance (see (16)). We obtain a similar result: 
a plateau for values of a around 0.8 is clearly visible. 
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FIG. 12. (Color online) Analysis of the stability with respect to the variation of the damping parameter from the value a — 0.85. 
The network is a scale free graph with 128 nodes generated using networkx for the quantum algorithm (colored lines) and of 
256 nodes for the classical one. a) In color the classical fidelity (see (12)) between the quantum PageRank calculated using 
a = 0.85 and the one calculated using values in the range from 0.01 to 0.98. In grey the one for the classical algorithm. 
We notice that the quantum PageRank is more robust with respect to the variation of the damping parameter. b)The trace 
distance between the quantum PageRank calculated using a — 0.85 and the one calculated using values in the range from 0.01 
to 0.98 (using (16)). Plot of the trace distance to compare the quantum PageRank obtained with a = 0.85 to the one obtained 
using other values of the damping parameter. 



the minimum fidelity between two values of a and a' does 
not go below the value of 0.91. One should compare with 
the analysis of the classical PageRank [51 , where the 
fidelity between different values of a and a' can be ap- 
proximately 0. It can thus be inferred that the ranking 
is rather robust when it is performed with the quantum 
PageRank. 

We have also investigated the behavior around the value 
of the damping parameter a = 0.85. From the classical 
fidelity between the quantum PageRanks at a = 0.85 and 
at a ranging from 0.01 to 0.98 one can see that there is 
a plateau around the value of a = 0.85 (see figure 11) 
extending especially for smaller values of a. There is a 
dip for a = 0.95 which is due to the fact that the ranking 
is very sensitive to changes in the damping parameter a 



when this approaches 1. For this value only the second 
term in the Google matrix G, giving random hopping, 
survives. 

The analysis was made more precise and the conclusions 
more cogent by considering also the measure of distance 
of rankings originating from the trace distance (see (16)). 



Also in this case the overall robustness of the ranking 
performed with the quantum PageRank is evident. One 



can see from figure 11 that the maximum value of this 
measure of distance is 0.18 comparing any two values of a 
and a' ranging in the aforementioned interval. It can be 
observed also that the region where this ranking is more 
robust is where a ^ 0.8. Indeed from figure pT] the blue 



region is wider (and correspondingly the ranking more 
robust with respect to perturbation of the value of a). 
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FIG. 13. (Color online) The fidelity obtained by applying 
the classical fidelity (see eq. (12)) in the case of the classical 
Page Rank. The network analyzed is a scale free network with 
256 nodes (generated with NetworkX). The damping param- 
eter varies ranging from 0.01 to 0.98. One can see that for 
extremal values of a, a' the value of the fidelity drops below 
0.4. 



We find also in this case that figure 12 the trace distance 
is rather smooth for a = 0.8. One can see that also in 
this case there is a curious peak for for a ~ 0.95 similarly 
to the previous case. 

We conclude that the quantum PageRank as measured 
by the classical fidelity or by the trace distance seems to 
vary mildly with respect to the variation of the damping 
parameter a. Indeed the minimum fidelity between any 
two distributions of importance arising from the quan- 
tum PageRank is 90%. The maximum of the distance 
obtained by using the trace distance between any two 
states (with different a) is 0.18. This means that the 
quantum PageRank is very robust with respect to varia- 
tion of the parameter that controls the fraction of random 
hopping. It is much more robust than in the classical 
case. From our analysis of the stability of the classical 
PageRank we find that for extremal values of a, a' the 
value of the fidelity is less than 0.4 (see fig.|V|). In ref. [51] 
the minimum fidelity was found to be very close to 0% 
between PageRanks' rankings corresponding to extremal 
values of the damping parameter a. 



VI. POWER LAW BEHAVIOR FOR QUANTUM 
PAGERANK 



In section [n] we have found that the quantum PageR- 
ank is able to highlight the structure of secondary hubs. 

In addition to this finding we have shown in section IV 
that this strength of the quantum algorithm is associ- 
ated with the fact that the quantum walker on networks 
of scale-free type is in a localized phase, i.e. the nodes 
with a significant average quantum PageRank are a neg- 
ligible fraction of the nodes in the network. On the other 
side, for the classical algorithm, it has been shown [55] [56] 
that for real networks the nodes' classical PageRanks /j. 



sorted in descending order, behave following a power law. 
That is, the classical PageRanks decrease approximately 
like Ij ~ where f3 ~ 0.9. This is a clear witness of 
the fact that the algorithm is able to identify the hubs. 
Furthermore, the scaling exponent (3 measures the rel- 
ative importance given to the hubs with respect to the 
other nodes. 

A similar study is therefore important in the case of the 
quantum PageRank and will be carried out in this sec- 
tion. Indeed, it is desirable that the quantum PageRank 
display a power law scaling behavior. Indeed, such a be- 
havior is distinctive of the fact that the algorithm be able 
to uncover the scale free nature of the network. 
In order to point out the scaling behavior of the quantum 
PageRanks of the nodes it is clear from the conjectured 
form 



(17) 



that we can extract such behavior from the slope of 
the log-log plot of the quantum PageRanks versus the 
(sorted) index of the i as can be seen from log{Iq{Pi)) ^ 

-Pq log(z) . 

Here we perform the analysis on scale free networks of 256 
nodes. We calculate the classical and quantum PageR- 
anks of the nodes and after having sorted the nodes in 
descending order we analyze the log-log plot of the clas- 
sical and quantum PageRanks versus the nodes' index. 
Considering one instance of a graph in this class one can 
clearly see (cf. figure 14) that both rankings obtained 
using the classical and the quantum PageRank display a 
power law behavior. The plot displays three areas, cor- 
responding respectively to the hubs, to the intermediate 
part of the list and the low part of the ranking (which in 
the classical case is degenerate in importance). 

The fact that this behavior persists also in the case 
of the quantum PageRank is due to the effect, similar 
to what has been observed in [51] for the classical case 
(and in section IV for the quantum Pagerank that the 
walker is in a localized phase in the case of scale free 
networks). Consequently the most highly ranked nodes 
tend to concentrate nearly the totality of the importance. 
This can be clearly seen in figure [l4| area /, where the 
hubs' importances lie above the line. 
Furthermore, it is clear that the scaling coefficients are 
different in the quantum and classical case. We find that 
Pq < Pel the quantum PageRank has a smoother behav- 
ior, giving less relative importance to the nodes in the 
high part of the list. On the other hand it is also able to 
better rank in the low part of the list (where the classi- 
cal PageRank gives highly degenerate values) lifting the 
degeneracy, (cf. the area /// in figure [l4|. 

The quantum PageRank is therefore able to better dis- 
tinguish the nodes' importances in the lowest part of the 
list. This is because the difference between the impor- 
tance of the nodes in the higher and lower part of the 
list is lower. Furthermore, the power law behavior in- 
terpolates a greater portion of the data with respect to 
the classical case as can be noted by inspection (the area 
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// in figure [14] extends much more in the case of the 
quantum algorithm). 

To complete the analysis we consider an ensemble of scale 
free networks in order to display the ensemble's proper- 
ties rather than the particular instance's. The ensemble 
consists of 29 scale free networks. It is clear that the 
properties found in the instance of the graph in figure 14 
persist also when considering a mean property of the en- 
semble (see figure 14). We conclude that these properties 
are generic and not an artifact of considering only one in- 
stance of the class of scale-free networks. 
Finally, we consider the real- world network, a subgraph 
of the WWW obtained by exploring pages linking to 
4T 



www.epa.gov 

ing (see figure 



which we refer as EPA in the follow- 
15|). Also in this case we find a power law 

Pel- 



behavior with j3q = 0.30 < 0.45 



The plot displays 



three areas as mentioned above for the other cases and 
our conclusions are valid also for this real-world network. 



VII. SENSITIVITY OF THE QUANTUM 
PAGERANK ALGORITHM UNDER 
COORDINATED ATTACKS IN SCALE-FREE 
GRAPHS 

In this section, we aim to study of the sensitivity of 
the quantum PageRank algorithm under coordinated at- 
tacks. More precisely, we ask how much the list of quan- 
tum PageRanks of an TV-node graph changes as a whole 
if the n most important nodes (hubs) of a network fail 
(e.g. due to a hacker attack) and the quantum quantum 
PageRank algorithm is run on the remaining {N — n) 
- node graph. Motivated by the fact that the real-world 
internet belongs to the class of scale- free networks, we fo- 
cus in our study on the scalefree networks of mesoscopic 
size (graphs of 16 and 32 nodes). 

Operationally, we proceed in our numerical study as 
follows: 

(i) First, we determine for the initial A/'-node graph the 
quantum PageRank values and the ordered list of nodes 
according to the quantum PageRank algorithm. 

(ii) Next, we take out the most important node (main 
hub) from the graph by eliminating the node itself as 
well as all its in- and outgoing links from and to the 
other (A^ — 1) nodes. On this reduced graph (with cor- 
responding modified connectivity matrix C) we carry 
out the quantum PageRank algorithm to determine of 
the modified list of quantum PageRanks, with a possi- 
bly modified order list of — 1 nodes. We note that on 
this modified graph the algorithm differs quantitatively 
from the one run on the complete A/'-node network, as 
the Hilbert space, the initial state and the underlying 
coherent dynamics are different due to the modified con- 
nectivity matrix C and resulting modified Google matrix 
G' (see Sec. [IB). 

(iii) We compare the ordered list of {N — 1) nodes ac- 
cording to the quantum PageRank values with the origi- 
nal list for the A/'-node graph, where the most important 



node is taken out (see Fig. [T6pi). To quantify the overall 
difference between these two lists of N — 1 elements, we 
employ Kendall's coefficient [57]. This function returns 
for lists in which the order of all elements is the same 
(irrespective of the individual values associated to each 
element of the list), zero for lists whose order of elements 
is exactly the opposite, and values in between for lists 
where the order of elements partially differs. 

This procedure of steps (i) to (iii) can be iterated to 
take out subsequently the n most important nodes ac- 
cording to the initial quantum PageRank list of the A/'- 
node graph. The resulting quantum PageRank list of the 
reduced {N — n) - node graph is then compared to the 
initial list (not including the n most important nodes). 

To compare the sensitivity of the quantum PageRank 
algorithm with the classical PageRank algorithm, we per- 
form the same type of coordinated attacks in the classical 
scenario, i.e., we analyze how the ordered list of classi- 
cal PageRank values changes when the n most important 
(according to the classical PageRank protocol) nodes of 
the graph are taken out and the classical PageRank al- 
gorithm is run to determine the importance of nodes in 
the reduced {N — n) graphs. 

The results are shown in Fig. [TgJd and c. We find that 
when the most important nodes (hubs) are attacked and 
fail, and the PageRank is recalculated for the reduced 
graphs, the order of the importance of the remaining 
nodes changes with respect to the list of the initial com- 
plete graph both in the classical and in the quantum 
case. In the quantum PageRank algorithm, attacks on 
hubs have a stronger effect than for the classical algo- 
rithm. This behavior can be understood by the fact that 
whereas for the classical case there is a large degeneracy 
of importance values of nodes of low PageRank, quantum 
fluctuations partially lift this degeneracy - see discussion 
in Sec.[ll|and the insets in Fig.|2j Thus, when hubs of the 
network are attacked and fail, the order of less important 
nodes - whose importance values slightly differ, can truly 
change in the quantum case, whereas the degeneracy of 
a larger number of nodes persists in the classical case. In 
other works, the increased capability of the quantum al- 
gorithm to resolve more structural details of the directed 
graphs, comes at the cost of an increased sensitivity to 
structural changes of the network. 



VIII. CONCLUSIONS 

In this paper we have studied the behavior of the quan- 
tum PageRank algorithm, developed in Ref. [3 , when ap- 
plied to complex networks. We have found that the quan- 
tum algorithm is able to clearly distinguish the struc- 
ture of the underlying network. More specifically, the 
quantum PageRank's behavior is distinctive for the three 
classes of complex networks studied in this work: scale- 
free networks, graphs of Erdos-Renyi type and hierar- 
chical networks. In particular, we have observed that 
the quantum algorithm when applied on scale-free net- 
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FIG. 14. (Color online) a) The plo t of t he logarithm of the PageRanks (upper part) and quantum PageRanks (lower part), (after 
being reordered, see text in section VI ) versus the logarithm of the node's label. As a guide to t he eye we have superimposed two 
lines with slope equal to —0.9. One can clearly distinguish three zones (see text in section VI). b) The plot of the logarithm of 
the means over the ensemble of graphs in a class of scal e fre e networks of the PageRanks (upper part) and quantum PageRanks 
(lower part), (after being reordered, see text in section \VJ\ versus the logarithm of the node's label. As a guide to the eye we 
have superimposed two lines with slope equal to —0.9 in the classical case and —0.85 in the quantum case. 



works is able to highlight the structure of the secondary 
hubs and to resolve the degeneracy in importance of the 
low lying part of the list of rankings, which represents a 
typical shortcoming of the classical PageRank algorithm. 
Although best suited for scale-free graphs the quantum 
PageRank is also able to univocally uncover whether 
graphs lie in the Erdos-Renyi class. Applied to hierar- 
chical graphs the algorithm has the capability to better 
reveal the hierarchy of levels, of which the graph is com- 
posed, and to highlight the connectivity structure within 
every hierarchy layer better than its classical PageRank 
counterpart. 

Regarding the quantum PageRank algorithm as a di- 
rected quantum walk, we have studied the localization 
properties of the quantum walker in the quantum pro- 
tocol. By an analysis of the Inverse Participation Ra- 
tio (IPR), we have observed localization of the quantum 
walker in the case of the quantum PageRank applied to 
scale-free networks under standard conditions (damping 
parameter a = 0.85). This finding is consistent the abil- 
ity of the quantum algorithm to highlight hubs of the 
network. In contrast, for Erdos-Renyi graphs delocaliza- 
tion was found, which is in accordance with the absence 
of hubs for this class of networks. 

Furthermore, we have analyzed the robustness of the 
quantum algorithm with respect to variations of the 
damping parameter. We find a very high degree of ro- 
bustness as compared to the classical PageRank protocol, 
which indicates that the value of this parameter, whose 
choice is to some extent arbitrary, turns out to be not 
crucial for the quantum algorithm to work reliably. 
Furthermore, we have found that the distribution of im- 
portance values of quantum PageRanks over the nodes of 
scale-free networks follow - as has been previously found 
for the PageRank - also a power law behavior. However, 
the corresponding scaling exponent is for the quantum 



protocol smaller than in the classical case, indicating a 
smoother ranking of nodes. In contrast to the classical al- 
gorithm, in the quantum protocol the hubs of the graphs 
do not concentrate the whole importance and the algo- 
rithm lifts the degeneracy of the large set of nodes with 
low importance values. This increased ranking capabil- 
ity comes at the cost of being more sensitive to struc- 
tural changes to the network such as coordinate attacks 
on hubs. 

Remarkably, the described characteristics of the quan- 
tum PageRank even persist if the algorithm is applied to 
real- world networks. We have studied and successfully 
tested the performance of the algorithm by applying it 
to a real-world network, originating from the hyperlink 
structure of www.epa.org [41 , thereby showing that its 
intriguing properties are not restricted to an application 
of the algorithm to artificially, numerically grown net- 
works. 

The classical PageRank algorithm has been the subject 
of exact studies yielding analytical results, and other in- 
teresting studies [58H6Q) . It would be nice if the quantum 
algorithm can yield also exact analytical results and for 
this the class of hierarchical networks is a good candidate. 

A related subject of current study is the entanglement 
properties of quantum complex networks [6TH63] . A dif- 
ferent line of studies has pursued the application of the 
quantum adiabatic algorithm to the classical PageRank 
algorithm [64 . 

In future work it will be interesting to analyze in more 
detail the impact of random failures of nodes in large net- 
works of differing topology. From an algorithmic point 
of view, it is an interesting task to develop a dissipative 
version of this algorithm and to understand its perfor- 
mance and robustness properties in such scenario. Dis- 
sipation has already been considered as an element with 
respect to some aspect of the algorithm [65l |66] , but the 
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FIG. 15. (Color online) Plot in log- log scale of the quantum 
and classical PageRanks of the nodes of a piece of a real net- 
work (from EPA). The quantum and classical PageRanks are 
displayed after hav ing been sorted in descending order, see 
text in section |VI[ ) versus the logarithm of the node's label, 
a) The classical PageRanks. As a guide to the eye we have 
superimposed the fitted lines (in log-log scale) with equation 
= ccir'^ where /3 = 0.4545 and Cd = 0.0185. b) The aver- 
age quantum PageRanks. In this case we have superimposed 
the fitted line (in log- log scale) with equation = Cqi~^ 
where /3 = 0.3066 and Cd = 0.0095. 



development of a truly dissipative version in the spirit of 
dissipative quantum algorithms and computation [671 EH] 
remains an open question. Furthermore the growing field 
of complex quantum networks would benefit from a ver- 
sion of the algorithm that is able to rank nodes in the 
more general case where qubits are located at the nodes 
of the network. An important question in this scenario 
is whether an algorithm based on a multi-particle quan- 
tum walk [69, 70 is needed in this context, or if there 
exists for this task an efficiently simulatable algorithm 
that belonging to the computational complexity class P. 
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FIG. 16. (Color online) Numerical study of the sensitivity 
of the quantum and the classical PageRank algorithm un- 
der coordinated attacks in scale-free graphs, a) Conceptual 
schematics illustrating the comparison of the PageRank list 
corresponding to networks, where the most important node(s) 
are taken out, with the PageRank list corresponding to the 
complete graph with all nodes intact, by means of Kendall's 
coefficient K. b) Numerical results for graphs of 16 and 32 
nodes, where up to 5 of the most relevant nodes have been 
taken out. The data has been obtained by averaging over 100 
random scalefree networks, which have been generated with 
networkX. Error bars correspond to one standard deviation. 
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APPENDIX 



From eq. (15): 



{Q{Pua) = (Xn (Tr2 [pl\2t)MP])) ^ (A.18) 

where for bookkeeping purposes in the derivations that 
follow it has been made explicit to which spaces the den- 
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sity matrix refers to. Using the fact that trace preserving 
quantum operations are contractive. That is: 

D {pl\2t),pl',{2t)) > D {Tv^pl\2t),Tnpl^2t)) 

(A.19) 

and making use of another property, namely that the 
trace distance of two states is attainable as the maximum 
over all strong (i.e. PVM) measurement outcomes of the 
difference of the two states: 



D{p,a) = maxTr[M {p - a)] 
that, in our case, specializes to: 



(A.20) 



maxTr2 

M(2) 



M(2) {TTrpi\2t) - Trip^?(2t))l (A.21 



we obtain Finally it is useful to measure the average dis- 
tance as we are interested in the stability of the average 
quantum PageRank: 



( max Tt2 

\m(2) 



{D {TT,pi\2t),TTrp'J{2t))) 
M(2) {TT,pi\2t)-TT,p'J{2t)) 



(A.22) 



eq. (A.22) can be rewritten as: 



max 

M(2) 



{D (Tripi2(2t),Tripi?(2t))) 
TriTr2M^2Va'(2t) - TtiTt2M^^^ pi^,{2t) 



(A.23) 



which in our case is clearly stated as: 



(D {TT,pi\2t),TT,p'J{2t)))^ = 
max I (/,(P,, m, a)) - {Iq{P,, m, a')) \ (A.24) 



Where we can take the absolute value because the dis- 
tance between two states is a nonnegative number. 
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